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Abstract 

We show that particle production in proton-nucleus (pA) collisions in the 
Color Glass Condensate model can be related to Deep Inelastic Scatter¬ 
ing of leptons on protons/nuclei (DIS). The common building block is the 
quark antiquark (or gluon-gluon) dipole cross section which is present in 
both DIS and pA processes. This correspondence in a sense generalizes the 
standard leading twist approach to pA collisions based on collinear fac¬ 
torization and perturbative QCD, and allows one to express the pA cross 
sections in terms of a universal quantity (dipole cross section) which, in 
principle, can be measured in DIS or other processes. Therefore, using 
the parameterization of dipole cross section at HERA, one can calculate 
particle production cross sections in proton-nucleus collisions at high en¬ 
ergies. Alternatively, one could use proton-nucleus experiments to further 
constrain models of the dipole cross-section. We show that the McLerran- 
Venugopalan model predicts enhancement of cross sections at large p± 
(Cronin effect) and suppression of cross sections at low px. The cross 
over depends on rapidity and moves to higher p± as one goes to more 
forward rapidities. 

SPhT-T02/167 

1 Introduction 

High gluon density effects at high energy (small x) have been the subject of 
intense theoretical and experimental investigation [1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 
12, 13, 14, 15, 16, 17, 18, 19, 20]. At small a;, the number of gluons per unit area 
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and rapidity in the wave function of a high energy hadron or nucleus becomes 
large and the high energy hadron or nucleus can be described by a classical color 
field. Since most of the gluons in the wave function of a high energy hadron or 
nucleus are in a coherent state characterized by their typical momentum Qg, a 
high energy hadron or nucleus is dubbed a Color Glass Condensate. An effective 
action and renormalization group approach has been developed which enables 
one to calculate cross sections in a high gluon density environment. 

In [21, 22, 23, 24], for proton-nucleus collisions at high energies and in the 
forward rapidity region^ we proposed describing the proton by the QCD parton 
model and the nucleus by the Color Glass Condensate model. One can then 
relate particle production in proton-nucleus collisions to the scattering of a 
quark or gluon from the Color Glass Condensate [21, 22, 23, 24]. Here we show 
that one can relate particle production in proton-nucleus collisions to Deep 
Inelastic Scattering of electrons (or more precisely, virtual photons) on protons 
and nuclei. This is due to the fact that the description of the target as a Color 
Glass Condensate is universal and independent of the process considered. 

While the relation between particle production in proton-nucleus collisions 
and DIS structure functions has been known for a while in the case of dilepton^ 
(virtual photon) production [25, 26] (and references therein), here we show that 
this relation is more general in our framework and holds for production of any 
particle (which has a known fragmentation function) in proton-nucleus collisions 
at high energy. This is crucial in view of the upcoming proton (deuteron)-nucleus 
experiments at RHIC. One would like to have predictions from QCD with as 
few model dependent assumptions about nuclear effects as possible. 

In proton-proton collisions, standard calculations are based on collinear fac¬ 
torization theorems, proven at the leading twist level. The predictability of the 
theory is due to the fact that the building blocks (distribution or fragmentation 
functions) used in the cross section are universal (i.e. process independent) ob¬ 
jects which can be measured in a given process and used to predict the outcome 
of other processes. It is well known that higher twist effects break collinear 
factorization and one has to resort to other methods. In proton-nucleus col¬ 
lisions, in order to take nuclear modifications such as shadowing and Cronin 
effect into account, one commonly modifies the parton distributions to include 
nuclear shadowing and introduces some model for Cronin effect. The parame¬ 
ters of the model are then constrained by fits to the existing data and used to 
make predictions for other processes or higher energies. 

However, it is clear that this approach cannot be true in general since 
collinear factorization theorems break down due to higher twist (multiple scat¬ 
terings) effects which become significant in scattering off nuclei. Therefore there 
is no guarantee that the parameters of such models are universal and process 
independent and these models have little predictive power. Here we show that 
one can generalize the standard leading twist expressions in such a way that 
all higher twist effects are included and the universality of the building blocks 

^By forward rapidity region we mean, roughly, any rapidity between proton rapidity and 
mid-rapidity. 

^We would like to thank Y. Kovchegov for pointing this out to us. 
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(dipole cross sections) of the theory is preserved and the predictive power of the 
theory is restored. 


2 Reminder: DIS and the dipole cross-section 

2.1 DIS amplitude 

We first start by a brief reminder of how the dipole cross-section appears in Deep 
Inelastic Scattering. This calculation will in fact be useful in order to compare 
with the very similar problem of dilepton (i.e. virtual photon) production in 
pA collisions. The nucleus is treated using the Color Glass Condensate model, 
i.e. as a set of classical and randomly distributed color sources that generate 
a classical field. There are three relevant diagrams for the DIS amplitude in 
this model, which are represented in figure 1. In this figure, the black dot 



Figure 1: The relevant diagrams for A ^ qqX in the Color Glass Condensate 
model. The black dots denote the eikonal interaction between the quark or 
antiquark with the classical color field. 

denotes the eikonal all-orders interaction between a quark or antiquark with the 
classical field of the nucleus. More precisely, for a quark or antiquark line in 
which the incoming momentum is P and the outgoing momentum^ is Q, the 
(time-ordered) eikonal scattering amplitude reads 

reik(Q,P) = 27r5(g--p-)7-sign(p-) J l) , 

( 1 ) 

where we assume that the nucleus is moving close to the speed of light in the 
+z direction. U{x±) is a matrix living in the fundamental representation of 
SU{N,), given by 

I7(a;x) = Texp| dz~^pa{z~, (2) 

with in the fundamental representation, and where pa{z~, z±) is the density 
of color sources in the nucleus. How to average over these sources is explained 
in [27, 23]. 

3 We use the following notations throughout this paper: P denotes a 4-momentum, p 
denotes its spatial part, pj_ its transverse components, and its longitudinal components 
in light-cone variables. 
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We first compute the A qqX amplitude, and bring it to a form that will 
be easy to compare with the photon production amplitude. A straightforward 
use of the CGC rules (see [27, 28, 23, 24] for details and examples) gives for this 
amplitude the following sum of three terms^, each corresponding to a diagram 
in figure 1: 




xu{q) 


- K + w)7 


(P) 


{Q — Ky — w? + it 

+i J (C/(a:m) - 1) 

7“(/f - Z’ + m)7^ 


xu{q) 


+i 


dHi 


(2^)2 

xu{q) 


{P — Ky — m? + ie 

(yxi±(yx2j_ 


(p) 


(U{xi^) - 1 ) {uyx2±) - 1 ) 


7 + to) 7''(^ -/f -/. + m)7” 


2p-[{q^-l±y + m‘^-2q-k+] + 2q-[{qj_-k±-l±y + m^] 


Kp) ■ 
(3) 


In this equation Q, P and K are the momenta of the outgoing quark, outgoing 
antiquark and incoming photon respectively, Xi 2 _ and X 2 ± are the coordinates 
in impact parameter space of the quark and antiquark. In the third diagram, 
the momentum L is the momentum transferred between the nucleus and the 
quark line. Its — component is zero because the nucleus is moving at the speed 
of light in the +z direction, and its + component has already been integrated 
out by using the theorem of residues (we pick the pole in the upper half-plane, 
at = 5 + — [(qj_ ~ ^J-)^ + iny — ie]/2q~). 

At this point, it is useful to note that since P and Q are on-shell, we have 
the following two identities: 

u(g) = ^M(q)7"(m -k , 

vip) =p)'y~v{p) . ( 4 ) 

2p 

Inserting them respectively in the first and second term of Eq. (3), and in¬ 
troducing also a dummy variable l± via f d^Xi±d^l±/(27r)^ exp(il± ■ Xi 2 _) or 
f d^X2±d^l±/(27ry exp(il± ■ X 2 ±), we can combine three terms into one: 

M'^^^(klq,p) ='-J J d^Xi2_d^X2± 

“^The factor 27r5{k~ — p~ — q~), which expresses the fact that the problem is invariant 
under translations in has been excluded from the definition of the amplitude 
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Xgiii-a.i_Lgi(p_L+qi-fc-L-Ji)-=«=2X {Uixi^)u\x2±) - 1) 

xu{q)T'^{k^,k±\q-,p-,q^^ -l±)v{p) , (5) 


where we denote 


rf^{k^,k_i\q ,p ,g_L - Zj.) = 

^ _ 7~(^ - /^ + w)7^(^ - /C - /. + m)7~ _ 

~ + m‘^-2q-k+] + q-[{qj_-k±-lj^)‘^ + m^] 

This is our hnal expression for the DIS amplitude off the target color held. Note 
that T'' a priori depends on and However, in its list of arguments, we 
have anticipated the fact that the p^^ and that comes from the spinors drops 
out in physical quantities, because in the squared amplitude they always appear 
in combinations such as j~u{q)u{q)j~ = 2q~. In Eq. (5), E'' is a quantity 
which describes how the virtual photon splits into a qq pair, which longitudinal 
momenta q~,p~ and transverse momentum of the quark q^ — l±, while the 
factor U{xl±)U^X 2 ±) — 1 can be seen as the scattering amplitude of the dipole 
on color held of the nucleus. 


2.2 DIS Cross-section 

From Eq. (5), one can obtain the ^*A cross-section as 


DIS 


(fiq (Pp 1 

(2^)22(70 (27r)32po2A:-^’"^^^ ^ ^ 

X {M^^^^^{k\q,p)MZs{k\q,p))^t^{K)el{K) , 


(7) 


where (• • ■)^ denotes the average over the color sources and where e^{K) is the 
polarization vector of the photon. At this point, in order to exploit the fact that 
depends only on m± = qj_ — C, we can use this quantity as the integration 
variable in and we can integrate out trivially the transverse momenta 

and pj^ of the quark and antiquark in the hnal state, in order to obtain^: 


dz 


327rfc2 


z(l-z) 


d'^r± 


d'^m± d^m'. 


(27r)2 (27r) 


2 


ik)e:ik) 


J d^X^ Tre (^1 - Re + ^)u\X^ - ^)) J 

xTrd ((</-f m)rf'{k^,k±\q~,p~,m±)(^- {k^,k±\q~,p~ , (8) 


with 2 the momentum fraction z = q /k {p =(1 — z)k ), where Tic denotes 
the color trace and Tr^ the Dirac trace, and where we have introduced the dipole 

^We have made explicit the fact that it is the real part of the correlator (■ ■ ■)p which 
appears in the total cross-section, in order to emphasize the connection with the optical 
theorem. However, the real part has no effect since this correlator is purely real in the Color 
Glass Condensate model. Therefore, we drop it in the following formulas. 
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size r± = Xi± — X 2 ± and barycenter X± = (a^i^ + X2±)l‘2.. If we introduce the 
dipole cross-section® 




r_L, 


adipoie(rx) = ^ / cPX^ Tr, (1 - C/(Xj. + ^)U\X^ 


(9) 


and the square of the photon wave function: 


xT:id{{(i+'m)V^{k^,ki^\q-,p~, m_ l )(^- m)V''^{k^,ki_\q~,p~, m'^)) , 

( 10 ) 


we can write the following well known formula for the 7 * A cross-section: 


J dzJ(fr±\'i/{k^,k±\z,r±)fadipoie{r±). ( 11 ) 


So far, we have treated the scattering of the quark and antiquark off the nu¬ 
cleus at a purely classical level, and disregarded its energy dependence. This 
energy dependence is usually encoded in the dependence upon the variable 
X = —K‘^f2K ■ Pn where P„ is the 4-momentum of a nucleon inside the nu¬ 
cleus"^. In the CGC model, this energy dependence is expected to arise in the 
functional which is used in order to perform the average over the sources, when 
one solves the renormalization group equation that controls its x-dependence. 


2.3 Models for the dipole cross-section 

2.3.1 McLerran-Venugopalan model 

A first possibility to estimate the dipole cross-section is to use the McLerran- 
Venugopalan model, i.e. to use a Gaussian distribution of the sources Pa{z~, Zd_) 
that contribute to the matrix U{xd^) defined in Eq. (2). Assuming a distribution 
of the form 


W[p\ = exp 


dx d^x± 


Pa{x ,x±)p^{x ,a;_L) '[ 
2p'^{x-,Xd_) j 


( 12 ) 


® Strictly speaking, the dipole cross-section should depend on the longitudinal momentum 
fractions .2 and 1 — 2 of the quark and antiquark, because these parameters control the rapidity 
interval between the quark or antiquark and the target. However, the small x evolution that 
we will consider later aims at resumming leading powers of ln(l/a:), which means that we 
neglect ln(l/ 2 ) and ln(l/(l — z)) in front of ln(l/x). This is a reasonable approximation at 
high energy (small x) and if the photon wave function does not overweight values of 2 close 
to 0 or 1. 

^At leading twist, and in the frame where the DIS can be seen as a quark of the nucleus 
being struck by the virtual photon, x has the interpretation of the longitudinal momentum 
fraction of the struck quark. 
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with , a;x) a density of color sources per unit volume, one obtains (see the 

appendix A of [27]): 


f^dipole (^_L ) — TtR 


1 — exp 



d'^z±[Go{z±) - Go(z± 



(13) 


where ^ a^ f dz~iJ?{z~) (in order to arrive at this formula for the dipole 
cross-section, one has to assume that the target is homogeneous in the transverse 
plane) and where Go is the free propagator in two dimensions: 


Go(^±) 


d^k± ^ ] f ^ \ 

(27r)2 kl "4^ ’ 


(14) 


where A is some infrared cutoff related to the scale at which color neutrality 
occurs. Therefore, A is at least as large as the inverse hadron size, i.e. A^^^. 
In fact, it can be proven that in a saturated target, color neutrality occurs 
over transverse spatial scales as small as [19, 29]. At small dipole sizes, 
this integral behaves like ln(l/r_LA). Quantum corrections to the McLerran- 
Venugopalan model give an x dependence to the saturation scale Qs and hence 
to the dipole cross-section [9, 10, 17, 18, 19, 20, 30, 31]. 


2.3.2 Golec-Biernat-Wiisthoff model - without evolution 

In [32, 33, 34], Golec-Biernat and Wiisthoff developed a model of dipole cross- 
section that is loosely inspired by the previous form of the dipole cross-section: 

CTdipoie(a;,r_L) = cto(1 - exp(-ri/4i?o(x))) , (15) 

with i?Q(a;)[GeV“^] = (x/xo)^- The main simplification compared to Eq. (13) 
is that the r± dependence is taken to be strictly Gaussian, i.e. one neglects the 
slowly varying logarithm in the exponential®. The x dependence is reminiscent 
of the BFKL evolution of the gluon structure function at small x. Then, in 
order to determine the values of the three parameters ctq, A and xq, they fitted 
the DIS data at HERA for x < 0.01. The best fit was obtained with (Tq = 23mb, 
A = 0.29 and xq = 3.10“^. 


2.3.3 Bartels-Golec-Biernat-Kowalski model - with evolution 

The previous model for the dipole cross-section was subsequently improved in 
[35] in order to make the small dipole limit agree with leading twist perturbative 
QCD, including the DGLAP evolution of the gluon density. This new model 

®When we connect particle production in pA collisions to the dipole cross-section in section 
4 , we will show that this Gaussian form fails to reproduce the standard perturbative result 
for the px spectrum of the produced particles. In other words, while this logarithm is not 
crucial in order to obtain a reasonable fit of small-a; DIS data, it is crucial in order to get the 
correct spectrum shape in pA collisions. 
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improves the agreement with the high data (up to < 450GeV^). In this 
version the dipole-proton cross section is [35] (see also [36]): 


^dipole(^_Lj — ^0 


- 3 ^^- 


(16) 


where = /Xq -I- C/r^^ and all the parameters are determined from fits to the 
DIS data from HERA. Using Eq. (16) in (11) gives the cross section for the 
interaction of a virtual photon with a proton target described by the Color 
Glass Condensate. Since we are interested in scattering off a nucleus target, we 
can replace xG{x,^^) in Eq. (16) by AxG{x, The parameter uo must also 
be scaled according to the transverse area of the target. 

The goal of the rest of this paper is to show that once one knows the dipole 
cross-section, there are several other interesting processes for which one can 
make quantitative predictions, since they can be expressed in terms of the same 
dipole cross-section. For instance, the cross section for jet production in pA 
collisions is very closely related to the dipole cross-section on a nucleus, as we 
will see in section 4. 


3 Virtual photon production in pA collisions 

3.1 Photon production amplitude 

The first example we consider is the production of virtual photons (or, equiv¬ 
alently, lepton pairs) in pA collisions, since this is in fact the process which is 
the most closely related to DIS. We follow here the description of the proton 
assumed in [23, 24], where it was assumed that the proton can be described in 
terms of standard parton distributions. Therefore, the relevant subprocess for 
the production of a virtual photon in pA is qA 97 *AT, which is related to 
7 *A ^ qqX by a crossing symmetry. 

The main difference in the case of 7 * production is the fact that there are 
only two diagrams, represented in figure 2. Note that there could a priori be a 




Figure 2: The relevant diagrams for qA < 77 *AT in the Color Glass Condensate 
model. The black dots denote the eikonal interaction between the quark or 
antiquark with the classical color field. 

third diagram with eikonal interactions of the quark both before and after the 
photon emission. However, this diagrams is in fact strongly suppressed by large 
energy denominators [23], which simply indicates that by the time the photon 



has been emitted after a first scattering of the quark off the nucleus, the nucleus 
is far behind and any further scattering is impossible (in order words, the nucleus 
is so fast that the photon is emitted outside of the nucleus). This discrepancy 
between the number of diagrams contributing to both processes is what makes 
the correspondence between the two non-trivial®, and worth checking explicitly. 
In the CGC model, the qA qj*X amplitude is given by^° 


M'^.{p\q,k) 



qj_+k±-pj_)-x± 


{U{x±) - 1) 


xu{q) 


+ /^ + to )7 
(Q + KY — m? + ie 


7 {Y - H 

{P — KY — rrP + ie 


u{p) ,(17) 


where the two terms correspond to the two diagrams of figure 2 respectively. 
In this equation p, q and k are the momenta of the incoming quark, and of the 
outgoing quark and photon respectively, and x± is the transverse coordinate of 
the quark. Using now the first of Eqs. (4) as well as 

u{p) = + p)-f~u{p) , (18) 

2p 

and introducing again a dummy variable l±, we can after some algebra rewrite 
the 7 * production amplitude as: 

AI((.(p|q,fc) = I J J (Yxij_(Yx2± 

Xe*U-XiXgi(9i+fcr-p_L-U)-a.2_L (U{xi±) - U{X2±)) 

xu{q) r^(-fc=^, -fcxig", -p~,q_L -1±) w(p) > (19) 


where T'' is again the function defined by Eq. (6). We see that the only difference 
between the DIS and 7 *-production amplitudes, besides the obvious changes 
P —P and K —K, is an extra factor UYx 2 ±) under the integral. In other 
words, the two amplitudes can be related by crossing symmetry up to a unitary 
matrix, which violates a strict crossing symmetry at the amplitude level. This 
has been noted independently by S. Peigne in [38], where Drell-Yan and DIS 
are compared up to the terms involving three scatterings on the target. In the 
present case, this violation is in fact due to the eikonal approximation. Indeed, 
this approximation being valid for an asymptotically large relative momentum 
between the projectile and the target, it does not allow to reverse continuously 
the momentum of the projectile. 

®If we were working at fixed order 0{as) instead of resumming all the eikonal interactions, 
only the first two diagrams of figure 1 would contribute to DIS, and finding the relation 
between DIS and 7 * production would be trivial. Note that a similar correspondence was 
found in [28], although in a case where the photon kinematics was far less general. 

Again, we do not include the factor 2tt5{p~ — q~ — k~) in the definition of ■ 
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3.2 Cross-section 

The photon production cross-section is obtained from the amplitude by 
d^q d^k 1 ^ 

= 7—— _ ,.3^, -— Iirdip —q —k) 

(27r)22(jfo (27r)^2fco 2 p 

x^jM>;,{p\q,k)M';t{p\q,k))^e^{K)e:{K) , (20) 

where the factor l/27Vc is for the average over the color and spin of the incoming 
quark. At this point, it is useful to replace the integration variable l± by 
m± = < 7 j_ — lj_ both in the amplitude and in its complex conjugate, in order to 
exploit the momentum dependence of the function T^. Another simplihcation 
comes from the fact that any term in the amplitude squared depends non- 
trivially only on two out of the four transverse coordinates. Therefore, the 
two transverse coordinates that appear only in exponentials can be integrated 
out immediately, which gives two delta constraints on some combination of 
transverse momenta, that can be used in order to perform the integrals over 
m± and Following this procedure, one obtains the following form for the 
differential 7 * production cross-section^^: 

j I ^ 

X j £x^Ti,(^l-U{Xj^ + ^-)UUXj_ 

Mip^- k±,pj_- k±) + M{qj_,q^) 

-M{qj_,pj_- k_i) - M{p^_- kj_,q^) , 
where we use the shorthand 

= e^{K)el{K) Trd((g^-k -p“, m^) 

{-k^,-k±_\q~,-p~. ( 22 ) 

In fact, the four terms in Eq. (21) could have been obtained directly by squaring 
the two terms of Eq. (17)^^. Again, we see that the dipole cross-section appears 
naturally in this cross-section. Therefore, if one determines the dipole cross- 
section from DIS data, then it becomes possible to make quantitative predictions 
for the production of virtual photons in pA collisions. This possibility has 
already been exploited in [26] in order to study nuclear effects on the Drell-Yan 
process. 

also use the freedom to add a constant to the term in UU^, because doing this does 
not change the final result for the cross-section. 

^^The form of Eq. (19) for the 7 * production amplitude is in fact not optimal in order to 
calculate the cross-section because it contains some redundant variables. Its only advantage 
is to make obvious the similarities with the DIS amplitude. 
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One can in fact obtain a form very similar to Eq. (11) for the 7 * production 
cross-section integrated over the transverse momenta and fcx of the final 
particles: 


cr^. = J dz j d'^r±\(j){p ,Pj_\k''~,z,r±)\^ adipoie{r±) 


(23) 


where z is the longitudinal momentum fraction taken by the photon (fc = zp , 
q~ = {1 — z)p~) and where we denote 


c/>(p ,p^\k+,z,r±)\^ = 


1 


647rp?_z(l — z) J (27r)^ J (27r) 




M(pj_ - fc_L,p_L - k±) + M{q^,q^) 
-M(qj_,p_L - k±) - M{pj_ - k±,q_p) 


(24) 


4 Forward pion production in pA collisions 

4.1 Quark scattering amplitude 

In this section we show how one can relate the cross section for forward particle 
production in pA collisions to the dipole cross-section. To be specific, we will 
consider pions but this can be repeated for any particle for which there is a 
known fragmentation function. We will use the results of [21, 22] for scattering 
of an on-shell quark from the target described by the Color Glass Condensate. 




> ^ > 


Figure 3: The relevant diagram for qA qX in the Color Glass Condensate 
model. The black dots denote the eikonal interaction between the quark or 
antiquark with the classical color field. 


The amplitude for scattering of a quark from the target described by a 
classical field is given by the diagram represented in figure 3, and has been 
obtained in [21, 22]. Factoring out the obvious 2Tr6{p~ — q~), this amplitude 
reads (the same as Eq. (1)): 


MqA^qx{p\q) 


u{q) 





u{p) , 


(25) 


where p and q are the momenta of the incoming and outgoing quark respectively. 
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4.2 Scattering cross section 

From this amplitude, it is easy to obtain the cross section [21, 22]: 

daqA^qX = (2:^^^27r5(p- - 9“)^ {MqA^gxip\q)M*qA^qxiP\Q))p , 

(26) 

where the prefactor l/2iVc comes from the average over the spin and color of 
the incoming quark. This can be rewritten as: 

X I d^XxTr, ( (uiXA + ^) - l) {uHX± - - l) -(27) 

Note that for the scattering of a gluon, one would have to replace the matrix U 
by its analogue in the adjoint representation, but the average over the color of 
the initial state would require a factor — 1) instead of l/Nc- Expanding 

the correlator in the bracket (• • •)p, we can rewrite this cross-section in terms of 
the dipole cross-section: 

daqA^qX = dq~6{q~ - p~) j 

X ^ J d^X^Ti,(2-UiXj_ + '^)-U\Xj_-'Y)) -^dipoie(r±) , 

(28) 

where, anticipating the use of collinear factorization for the incoming quark, we 
have set pj_ = 0. The first term (2 — U — is essential in order to obtain 
the correct total cross-section. Indeed, integrating over qj_ generates a d(r_L), 
and since crdipoie(’'± = 0) = 0, we have for the total qA qX cross-section^^: 

= ^J d^X^Tr, (2 - U{Xj_) - U^X^))^ . (30) 

If we recall (see [27] for instance) that {U{x±))p is suppressed exponentially 
like exp(—Qg/A^^^) when Qs A> A^^^, we can approximate the qA qX 
cross-section by: 

72 p 

daqA^qX = dq~S{q~ -p~) J d?r- (Jdipoie(r'_L)) , 

_ (31) 

order for the optical theorem to hold, the forward elastic (qA qA) amplitude must 
be; 

= 2p- I (F(*x) - l>p . (29) 

In other words, the forward qA —> qA amplitude is obtained from Eq. (25) by setting p = q 
and by averaging the amplitude over the classical color sources in the nucleus. 
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for a target of radius R. The term in 2'kR^ contributes only to the scattering 
in the forward direction (i.e. at qj_ = 0). 


4.3 Discussion 

4.3.1 Cronin effect at the partonic level 


From the previous relations, it is trivial to write the non-forward part of the 
-spectrum as: 


—>qX 

(P(l± 


1 


^dipole (^_l) 5 


(32) 


where (Tdipoie(9_L) is the Fourier transform of the dipole cross-section, defined 
as: 


<^dipole(Q±) — 


d^r_Le*'*^ '^^CTdipoie(r±) . 


(33) 


In other words, the g_L-spectrum of particle production in pA collisions is given 
by the opposite of the Fourier transform of the dipole cross-section. 

In the McLerran-Venugopalan model, this Fourier transform has been stud¬ 
ied in detail in [27, 28] where it controls the production of qq pairs in ultra¬ 
peripheral AA collisions and was denoted 'KR?C{qj^) (see for instance the ap¬ 
pendix B of [27]). In particular, the large q±^ behavior of this function was found 
to be 


C{qi) 







(34) 


Note that if one replaces by its expression in terms of parton distributions 
inside the target, the first term in l/q\ reproduces the perturbative QCD result 
for qg qg scattering with one gluon exchange in the t-channel (i.e. the 

dominant piece at high energy). In addition, if one assumes that scales 
like 4^/^ for large nuclei, then the leading term of the q± spectrum simply 
scales like A (another A?!^ comes from the transverse area ttR?). If the second 
term in this expansion scales differently with A, then there is a Cronin effect. 
From Eq. (34), one can see that the next-to-leading term in the g_L-spectrum is 
scaling like 4^/^ and is positive. Therefore, one has a positive Cronin effect in 
the McLerran-Venugopalan model. The ratio {A-ceilA)duqA^qXIdoqA^^t^qX, as 
predicted in the McLerran-Venugopalan model, is plotted as a function of q±^ in 
figure 4. The reference nucleus is chosen such that the corresponding saturation 
scale is Q^ref = 1 GeV^, and the ratio is displayed for two nuclei such that the 
saturation scales are Ql = 2 GeV^ and = 3 GeV^. One can see that this ratio 
goes to 1 at large transverse momentum, and exhibits a pronounced maximum 
at intermediate q±A before dropping below 1 at small transverse momenta. 

Alternatively, one could take Aref = A = 200 and consider quark nucleus 
scattering at mid rapidity, assuming that Q^ref = 1 GeV^ at mid rapidity. Then 
the solid and dashed lines correspond to the relative enhancement of the cross 
section in the forward rapidity region, y ~ 2.5 and j/ ~ 3.5 units away from 
mid rapidity (towards the proton), respectively. The low q± spectrum is more 
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= 1 GeV^ (A = 200 MeV) 



Figure 4: Cronin effect at the partonic level in the McLerran-Venugopalan 
model. The reference nucleus is chosen such that the corresponding satura¬ 
tion scale is Qlrei = IGeV^. This reference value of Qs is compared to nuclei 
for which Qg = 2 GeV^ and = 3 GeV^. The infrared cutoff is set to 
A = 200 MeV. 


suppressed as one goes to more forward rapidities and high q± enhancement is 
stronger and moves to higher q±. 

It would be interesting to study whether this phenomenon is present in 
more sophisticated models of the dipole cross-section, like the model of Bartels, 
Golec-Biernat and Kowalski. It is also necessary to fold this “partonic level” 
calculation with the proton parton distribution functions, and with the appro¬ 
priate fragmentation functions, as this could potentially suppress the effect. A 
numerical study of these issues is underway. 

4.3.2 Further constraints of the dipole cross-section in pA 

In DIS, one does not probe thoroughly the dipole cross-section, but only its 
overlap with the square of the photon wave function. This implies that some 
values of r± matter more than others in the integral. One could therefore see 
Eq. (32) as another way of constraining the dipole cross-section. In particular, 
one would like to recover the results of perturbative QCD at large q±, that is a 
cross-section that falls like qj^ (up to logarithms coming from the running Og 
and from the DGLAP evolution of the gluon distribution inside the target). 

For instance, this remark pretty much excludes the Golec-Biernat-Wiisthoff 
model of the dipole cross-section. Indeed, in this model (see Eq. (15)) the Fourier 
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Figure 5: Values of —(7x'^dipoie(9±) as a function of q±, for several values of x 
in the model of Bartels, Golec-Biernat and Kowalski. 


transform that gives the ( 7 j_-spectrum can be performed analytically, giving: 

dCqA^qX ^ Ro{x)cro 
cPqx^ TT 

The major difference compared to the result in the McLerran-Venugopalan 
model is the large qx behavior, which exhibits a Gaussian tail, as opposed to 
a power law tail. This fact alone is probably enough to justify not considering 
this model any further in order to study the region of large transverse momenta 
in pA collisions. 

For the Bartels-Golec-Biernat-Kowalski model, one needs to compute nu¬ 
merically the Fourier transform of the dipole cross-section. In order to do this, 
we have taken the parameters of their second fit (“Fit 2” in Table 1 of [35]), 
which lead to the results displayed in figure 5. One can see that for very small 
values of a; (x = lO”"*, x = 10“® on this plot) one obtains the expected scaling 
at large qx- Indeed, q\_ times the Fourier transform of the dipole cross section 
is almost constant. However, for larger values of x (x = 10“^, x = 10“^ on the 
plot), the Fourier transform of the dipole cross-section suddenly drops at some 
q± and its sign changes 

Above the qx where this happens, one would obtain a negative qA —> qX 
cross-section from formula Eq. (32). In this formula, one apparently probes 
different aspects of the dipole cross-section compared to DIS, and the model by 
Bartels, Golec-Biernat and Kowalski, although very successful at fitting all the 

^^This depends on whether a running or constant Os is used [39]. For constant Os one 
always gets a positive sign. See also eq. 20 in [35] and the discussion afterwards. 
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HERA DIS data at a; < 10“^, does not give a consistent answer for pA collisions 
at moderate x. Somehow, this defect is hidden in DIS by the way the photon 
wave functions weights the different values of Given this, one could perhaps 
reverse the main argument of this paper, and consider pA collisions as a way of 
either fine tuning the dipole model, or as a way of ruling it out. 


4.4 Prom qA qX to pA tiX 

To relate (31) to proton-nucleus collisions, we can use collinear factorization on 
the proton side and convolute (31) with the quark distribution function inside 
a proton and a quark-pion fragmentation function to get the cross section for 
pA 'kX. 


dcTpA — ^7T{k)X 

dk-cPkx 


J dXq dzqp{xq) 


dcTqA^qX 

dq~d?qA 


Dq/Tr{z) 


(36) 


where qp{xq) is the distribution function of quarks with fractional momentum 
Xq inside a proton and DqjAA) is the fragmentation function of a quark into a 
pion carrying fraction z of its energy. Using (31) in (36) we get 


dcr. 


pA^-iT{y,k_i_)X 




dyd^kx (27r)^ 

X J - crdipoie(r_L,a;3)) 


dzqp{xq) Dq/Az) 


(37) 


where we have used the following kinematical relations Xq = jZ\fs^ Xg = 

k±e~'^jZy/s, z„iin = kxt^l\fs [a/s is the center of mass energy for a proton- 
nucleon subsystem)- This way one can relate production of pions, kaons, pro¬ 
tons, neutrons, photons and dileptons in proton-nucleus collisions in the forward 
rapidity region to the structure functions measured in DIS. For this approach 
to be self-consistent, one needs to make sure that the dominant contribution 
to (37) comes from the small x (< 0.01) region. Our preliminary result seem 
to indicate that this is indeed the case. A detailed numerical study of particle 
production in proton-nucleus collisions at RHIC is currently underway and will 
be reported elsewhere. 


5 Conclusions 

In this paper, we have shown that in a classical description of the gluon con¬ 
tent of the nucleus, there are many interesting physical quantities that can be 
related to the dipole cross-section (or, at a more formal level, to the correlator 
(^U{0)U^{xx)) of two Wilson lines). This was already known for deep inelastic 
scattering and for the Drell-Yan process when one considers cross-sections inte¬ 
grated over the phase-space of the final state. In this paper, we establish also 
a similar correspondence for forward particle production in proton nucleus col¬ 
lisions, this time for the differential cross-section. More precisely, we show that 
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the q± spectrum of the produced particles is related to the Fourier transform of 
the dipole cross-section. 

At a more formal level, we have generalized the standard leading twist 
collinear factorization approach to calculation of hadronic cross sections by 
showing that there is a universal object, the quark antiquark dipole cross section, 
which appears in all particle production cross sections in high energy proton- 
nucleus collisions. This dipole cross section plays the role of leading twist parton 
distributions in an all twist environment. In principle, it can be measured in 
DIS or for example, in single inclusive pion production in high energy proton- 
nucleus collisions. Once it is measured in a given process, it can be used to 
predict cross sections for other processes such as single inclusive hadron, jet, 
photon or dilepton production in high energy proton-nucleus collisions. 

This result could therefore be used in order to make predictions for pA 
collisions, or to further constrain (or to rule out) the dipole model with proton- 
nucleus experiments. The detailed numerical study needed in order to make 
quantitative statements is in progress and will be presented in a future work 
[40]. However, a preliminary investigation of the Fourier transform of the dipole 
cross-section in the model of [35] seem to lead to a negative qA qX differential 
cross-section for some values of x at large qj_. This suggests that the small r± 
improvement of the dipole cross-section made in [35] over [32, 33] is still not 
quite able to ensure that the results of perturbative QCD are recovered at large 
transverse momentum. 
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